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1. Introduction 

The theory of holonomies and geometric phases associated with evolutions of a quantum 
system is by now a well developed subject. The initial work by Berry [1] on the 
Abelian geometric phase of adiabatic evolutions of non- degenerate states has been 
extended in many directions, to include non-adiabatic evolutions |2J as well as non- 
Abelian holonomies of sets of degenerate states [31 0] , and to mixed states [5] . It was 
subsequently discovered that the quantum geometric phase had an early counterpart 
in the geometric phase discovered by Pancharatnam in the context of classical optics 
[SI |7J [8] . The Pancharatnam construction has been generalized to the non-Abelian case 
by utilizing subspaces [9j [10J [TTJ [12]. More recently holonomies that bear a relation to 
correlations have been constructed in the context of multipartite and lattice systems 

[131 HH HH1 HE]. 

The idea of this paper is to develop the concept of geometric phases in another 
direction and use two-particle interferometry to construct correlation induced non- 
Abelian holonomies in a way that does not depend on degeneracy, but on the ability to 
divide the system into spatially separated subsystems. Instead of considering parallel 
transport of subspaces, we consider the natural tensor product structure of a bipartite 
state, induced by the spatial separation of the two subsystems, and local unitary 
operations to define the parallel transport. 

The parallel transport condition to be introduced is similar to that of the 
Pancharatnam construction [HI d [HJ , in which two states | A) and | B) of a quantum 
system are defined to be in-phase if their scalar product (B | A) is a positive number. 
This condition can be implemented in a Mach-Zehnder interferometric setup, where the 
spatial state of the system prior to the last beam splitter is a coherent superposition of 
the two paths [TTJ [18]. If we let | A) and | B) be the internal states corresponding to 
the output of respective paths of the interferometer, then Pancharatnam parallelity is 
achieved by shifting a U(l) phase in one of the paths, such that the interference intensity 
is maximal. 

In two-particle interferometry [191 1201 HI] the spatial state of the two-particle 
system is a coherent superposition of two distinct pairs of correlated paths for the two 
subsystems. In the same spirit as Pancharatnam we use a two-particle interferometric 
intensity, namely the coincidence intensity in a Franson interferometer [201 [ZD [22], to 
define an "in-phase" condition and the corresponding parallel transport. An arbitrary 
unitary operation is performed in one of the arms of the interferometer on the first 
subsystem, thus making the outputs of the two possible pairs of paths different. 
Subsequently, another unitary operation is performed in the other arm on the second 
subsystem, and is chosen such as to achieve maximal coincidence intensity. This second 
unitary is considered to be the "phase" degree of freedom of the system and at maximal 
intensity the two outputs are considered to be "in-phase", or "parallel". Thus we 
consider the orbit space formed by the state space of the system modulo this unitary 
degree of freedom on the second subsystem to be the space in which the system is 
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parallel transported. 

In the special case of pure two-qubit states, and evolutions generated by local 
SU (2) operations, the state space naturally fibrates through the second Hopf fibration 
and the orbit space of the two qubit-states can be mapped to the state space of a 
quaternionic qubit [23]. The coincidence intensity of the Franson interferometer, in 
this case, corresponds to the quaternionic quantum mechanics analogue of the Mach- 
Zehnder intensity and the associated parallel transport condition corresponds to the 
Levay connection [21] restricted to local SU(2) evolutions. 

The outline of the paper is as follows. It turns out that the Stokes tensor formalism 
[25] is convenient for our analysis. Therefore we briefly review this representation 
in section [2j Section [3] contains a description of the Franson interferometric setup 
and we define the parallelity condition in this setting. In section IH we describe the 
parallel transport procedure, discuss the properties of the related holonomy group and 
introduce the corresponding connection form. Finally, in section [5] we consider the 
parallel transport scheme in the special case of pure two-qubit states, and evolutions 
generated by SU(2) operations and its relation to the quaternionic representation of 
pure two-qubit states and the corresponding Levay connection. The paper ends with 
the conclusions. 

2. Stokes tensor formalism 

In the Stokes tensor formalism [25], single particle quantum states are represented as 
real vectors and iV-partite states are represented as real iV-tensors. A multi-partite 
system consisting of parts A, B, . . . Z, where part K has dimension D^, is represented 
as a D\ x D\ x . . . x D\- dimensional tensor Sj A j B ,,,j z . This tensor is related to the 
density matrix representation p of the same state as 

s 3A]B .., z xf A ®xf B ®---®x^ 

(1) 

where Xo = 1 X and xf K (jK = 1)2, ■■■,Dj ( — 1) are the D\ — 1 traceless generators 
of U(Dk) operations on system K. In the following, we shall in most cases use 
the simplifying notation xf f° r the generators on subsystem K. The Hermitean, 

traceless and linearly independent generators {xf}j=i of U(Dk) satisfy orthogonality 

Tr(xfxf) = 25 jk , and [ X f,xf] = 2i Zfii' fkimXl as well as {xf,xf } = ^5 M 1 K + 
d' 2 —1 

2 Xln=i dkinXn j "where and fju are the symmetric and antisymmetric structure 
constants, respectively, of U{D K ). The factors [o~oj K {D K — 2) + 2] _1 are inserted so 
that S jAjB ... jz = Tr (pxf A ®xf B ®- ■ -®xfj- 

We shall represent unitary operators U(Dk) on a form compatible with 
this formalism. These are represented by complex vectors with elements Uj = 



=££•••£ 

.74=0 7b=0 7z=0 



n 

K=A 



5 0jK (D K -2)+2 
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[MDk-2)+2] Tr (^f ) S ° that 

D 2 -1 

# = E Utf. (2) 
i=o 

Unitarity of the operators demand that {Uj} are complex numbers satisfying 

W+TT E l^'| 2=1 ( 3 ) 



and 



_D 2 -1 



j,k=0 

for each 1 < I < D 2 K - 1. 

In this paper we focus on bipartite systems A + B and therefore it can be instructive 
to consider some general properties of the Stokes 2-tensor. The zeroth row and zeroth 
column of the tensor, with elements Soj and Sjo respectively, are the Stokes 1-tensors 
corresponding to the reduced states of subsystem A and B, and these contain all the 
local information of the bipartite system. The remaining part, formed by the elements 
Sij, i,j > 0, contains all the correlations of the state and this subtensor, the correlation 
matrix, is denoted by M^. This matrix has previously been used to study correlations 
and separability in quantum systems [261 [27] . 

A unitary transformation on subsystem A transforms the reduced density 
matrix as Up A U^. This corresponds to a transformation RS, on S where Rjk = 
((d a -2)s 0+2) ^i&Xk U^xf )• The R matrix is an orthogonal matrix, as can be seen 
by observing that 

Dl-l D\-l 

[{D A -2)6 to +2\[{D A -2)6 j0 +2\j-£ ^ 

(5) 

Similarly a unitary transformation on subsystem B corresponds to an orthogonal matrix 
acting on S from the right. It should be noted that the M matrix transforms under local 
unitary operations on subsystem A and B through left and right action, respectively, 
by orthogonal matrices. 

As an example of how the M matrix behaves we can consider a pure two-qubit 
state in the Schmidt basis 

|^) = q 1 00) +6|11) (6) 

where a and b are real non-negative numbers and a 2 + b 2 = 1. Expressed in the 
Schmidt basis, the correlation matrix with elements Mjk = Tr ypaf®^), o~f and a^, 
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j, k = 1,2, 3, being the standard Pauli operators on subsystem A and B, respectively, 
reads 





( c 








\ 


M = 





-c 









V 





1 


/ 



where C = 2ab is the pure state concurrence [28]. Here it is clearly seen that M is 
rank one for product states. It must be emphasized that Mn = — M22 = C because 
we choose a, b G K. The elements of M are not explicit functions of concurrence 
for arbitrary complex coefficients a, b. However, since all pure two-qubit states with 
the same concurrence can be related by local unitaries, corresponding to orthogonal 
transformations acting from the right and the left on M, the absolute value of the 
determinant of M 

|detM|=C 2 (8) 

is invariant under local unitary transformations and measures concurrence. For 
maximally entangled states we thus have that |det M\= 1. 

When we consider mixed states there can be correlations also in separable states. As 
an example of how the M matrix registers correlation for mixed states we can consider 
the Werner states (29] 

Pw = pm^\+ { -^-i A ®i B , (9) 

where \ip) is some maximally entangled state and p G [0, 1]. The absolute value of the 
determinant is | detM |= p 3 , which can be compared to the square of the concurrence 
. Thus for p < I the determinant of M is nonzero even though the 
state is separable. This underscores that the M matrix for mixed states is sensitive to 
correlations in general. 



C 



max 



(3P-1) 



3. Parallel transport condition 

3.1. Inter jerometric setup and parallelity 

A Franson interferometer [2T)l [2T], I2"2"].is a two-particle device composed of two identical 
unbalanced two-path interferometers, so-called Franson loops, as shown in Fig. 1. The 
two parts of the bipartite state are emitted, one into each Franson loop. A beam splitter 
divides each path in two different paths of unequal length, which later converge at a 
second beam splitter. The difference in path length between the two arms in each 
Franson loop is chosen to be the same, and such that the difference in transit time At 
is greater than the single particle coherence time. The transit time difference between 
the paths must be smaller than the coherence time of the bipartite state to allow two- 
particle interference. Of importance is that the emitter is such that is impossible to 
define a time of emission. Detectors are placed after the convergence of the two paths 
in each Franson loop and coincidence measurements are made. 
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V-D 



source 



Figure 1. The Franson interferometer setup. Unitary operations on the internal 
degrees of freedom the particle are performed in the long arms of each Franson loop. 



Since we discard non-coincidental detections, corresponding to the bipartite system 
traversing a long path in one of the Franson loops and a short path in the other, it is 
necessary that the time resolution of the detectors is smaller than At. Furthermore, due 
to the requirement that time of emission cannot be defined, and since no measurements 
are made inside the Franson loops, it cannot be ascribed to a coincidence detection 
event that the bipartite system traversed either the two short paths or the two long 
ones. Hence the system is in a coherent superposition of having traversed the two 
long arms, and having traversed the two short arms. In the long paths of each sub- 
interferometer we place devices that perform unitary operations on the internal state 
of the bipartite state. After the point of convergence of the two paths, the effective 
unnormalized internal state is under the above requirements therefore 



p= _(i + U®V)po(l + 
where po is the initial internal state. Given U ■ 



(10) 



E J= o tyx 



lA f, and V 



B 



where xf an d xf are the generators of U(D^) and U(D B ), respectively, the coincidence 
detection intensity I AB is 



rAB 



The coincidence detection intensity I AB , henceforth referred to simply as "intensity", 
is the ratio between the measured intensity for a given U and V, and the intensity 
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measured if U = V = 1. The expression for I AB in the Stokes tensor formalism is 

D 2 _ l D 2 _ l 

iAB = \ + \Y. E MViUjSn (12) 

k=0 j=0 

where the second term is the interference term. 

We now define the parallelity condition in the Franson setup for a bipartite system 
consisting of two qudits of dimension Da and Db- We ask, given that a specific 
unitary operation U G U(Da) has been chosen in the first Franson loop, what unitary 
operation V G U(Db) should be chosen in the second Franson loop in order to maximize 
the coincidence intensity? We take maximal coincidence intensity as the definition of 
parallelity between the output of the two short paths and the output of the two long 
paths. This maximization procedure is the analogue of the procedure used to define 
Pancharatnam parallelity in the context of a Mach-Zehnder interferometer [T7J[I8], but 
here the Franson coincidence intensity has taken the role of the Mach-Zehnder intensity 
and V has taken the role of the U(l) phase factor. It should be noted that if S does not 

D 2 — 1 

have full rank, then there exist U such that X]fc=o U k Skj = for all j. In this case the 

D 2 — 1 

interference term is identically zero for all V, but if X]fc=o U k Skj ^ there will always 
be a V corresponding to maximal intensity. 

To find a formal expression for the operator V that maximizes the intensity, we seek 
to maximize I AB in equation (fl~2l) with respect to the coefficients V k , using Lagrange's 

D 2 — 1 

method. To enforce unitarity of V we introduce the constraints \ V \ 2 +-jj^J2j=i I 



vD£-l 



for each 1 < I < D 2 B 



Vj\ 2 = 1 and Ei,fc=o V j V k [ d jki + ifjki + (Sjohi + hoSji)] 
We thus construct the auxiliary function f ({Vj} , {V*}) , that is to be extremized, as 

D\-1D%-1 



fc=0 
D 2 -l 



3=0 



-A| V V* + 



D, 



E m 



Dl-l Dl-1 



v i v k [djki + ifjki + (S j0 5 k i + 5 k0 Sji)) 

1=1 j t k=0 



(13) 



Using Lagrange's method we seek the points were the gradient of the auxiliary function 
with respect to the variables Vj and V* vanishes. The components defining these 
points satisfy the equations 



A 



D 



B 



1 



D 



B 



>k() 



i=i 



D\-l 



(5 H V + S k0 Vi) + ^ v i( d jM + ifjki) 



3=1 



E U j S 3k, 



0<k<D 2 B 



1. 



(14) 
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If the coefficients Uj are ordered as a Z^-dimensional vector u and likewise the 
coefficients Vk are ordered as a D^-dimensional vector v, the above equations can be 
reexpressed matrix equation 

Bv = -S T u\ (15) 

where B is a A and dependent Hermitean matrix, given by 



By = X 



Z0°fc0 



+ fijidjkfilO + <5jZ^0fc + — ifjkl)- 

i=i 

(16) 

Provided B(X,jl) is invertible, the formal solution for V" can be given as 



D|-l 



v-= B- k 1 (x,fi)s^u;xF, (i7) 

l,k,j=0 

where ft = {//■/}. The explicit form of the Lagrange parameters, and thus -B -1 , is 
found by solving for the unitarity constraints on V. The solutions of these constraint 
equations give us the critical points of the intensity as a function of U. We can see from 
the constraints that for each solution A, \x there is a solution —A, —p, and if one of them 
corresponds to a local maximum, the other corresponds to a local minimum. There 
will be a unique solution to the maximization problem if and only if there is a unique 
global maximum of the intensity as a function of U, corresponding to a combination of 
parameters A,/2 satisfying the constraints. 

In the general case, finding this solution as a function of U appears to be a non- 
trivial problem. For product states, however, the unitary V that maximizes the intensity 
is easily found and is always a Abelian U(l) phase factor. This can be seen by observing 
that for product states p A ®p B the coincidence intensity is 

I AB = \ + ^ReTr(f)p A ) Ti(Vp B ). (18) 

and therefore the V that maximizes this expression is found to be V = e~* arg ( Tr ( f//3A ^i. 
We may note that for the trivial case where U = e^l, we see that the intensity is 
maximal if and only if V — e~^i, regardless of the state of the bipartite system. 

3.2. Example I: Qudit-qubit 

In section 13.11 the solution V G U(Db) to the maximization problem is not given 
on a closed form and it is not apparent how to find it. However, for the case 
where the second subsystem is a qubit, and therefore V G U(2), the B matrix is 
By = XSik + Y^j=i^i [{fijkfiio + <5j2<5fco) — ifjki]- If can be seen that B now separates as 
B = Xl + H, where H is Hermitean and H 2 oc 1, and therefore B^ 1 = -^==(— Al + H). 
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Explicitly 



B- 



Vdetfi 



/ -A 
I'l 

/'2 
\ /^3 



I'l 

-A 

-Z>2 



/'-2 

-i/j, 3 
-A 



A*3 

-A 



det£ = (-A 2 + /2-/i)' 



(19) 



While it is still not obvious how to find the Lagrange parameters in general we may 
solve it in some special cases. To illustrate this, let us consider a U that has the form 
U = U\0\ + U 2 a 2 , where <j\ and a 2 are the standard Pauli operators, and a pure state 

| V> = a 1 00> H- 6 J 11), (20) 

where a, b > 0. As a consequence of our choice of basis S w = S 20 = and the 
concurrence C = 2ab. This together with our special choice of U, implies that the 
intensity is 

jAB 



~ + jMViUi - v 2 u 2 ) 



(21) 



Using equations ( fTT|) and ( jT9j) . and the constraints, we find that p, = 0, A = ±C, where 
the positive sign corresponds to the global maximum while the negative sign corresponds 
to the global minimum. The unitary operator V corresponding to the maximum is 

V = U*a t - U;a 2 , (22) 

and the maximal intensity is 

1 



rAB 



;i+c). 



(23) 



3.3. Example II: Restriction to SU(D) 

A variation of the maximization procedure is to restrict the set from which U and 
V can be chosen. One natural restriction would be to consider only SU(Da) and 
SU(Db) operations in the Franson loops. When this restriction is made the qualitative 
properties of the parallel transport may change. It is for example no longer obvious that 
the unitaries V associated to a product state, will be a commuting set in the general 
case. The restriction where the second subsystem is a qubit and therefore V G SU(2), 
however, leads to a significant simplification of the maximization problem. Here, we 
solve this problem and in particular show that product states are indeed associated 
with commuting sets of unitaries. 

Since SU{2) can be parametrized by four real numbers subject to only one 
constraint, the solution of the maximization problem can be found easily for arbitrary 
states and arbitrary U e SU(Da)- We choose the parametrization of V such that Vo is 
real and i\ r i 1 iV 2 ,iV 3 are purely imaginary. The intensity in this parametrization is 



rAB 



1 1 

2 + 4 



Di 



=0 k=l 



(24) 
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and V is found to be 



v = Tx 



D\-l 



D\~l 3 



E & + u D s >oi B + i J2 E(^- - *w 

j=0 j=0 k=l 



(25) 



where cx s are the Pauli operators. The remaining Lagrange parameter A is found from 
the unitarity condition ^ _ I l 2= ^ an( ^ xs 



A = ±- 



U 



Di-l 



E + *w 



./0 



fc=i 



D 2 .-l 



.A- 



i=o 



(26) 



The sign of A must be chosen positive since the trivial case U — 1 implies V" = 1. 



When the bipartite state is a product state p 



P J 



® p B we find that V 



i[E?= A o Vj + U^)S jQ \ B + 2^ Efio ^ - U*)S j0 (p B - I s )]. Hence for product states 
the unitary V that maximizes intensity commutes with p B for any U. Therefore the 
set of unitaries V associated to a product state commute with the density operator and 
with each other. The signature of a product state is thus that the corresponding set of 
unitaries will be commuting when the unitary operators on the second subsystem are 
restricted to 577(2). 



3.4- Example III: SO(D) and two-redit states 

Another variation of of our procedure is to maximize the intensity for V G SO(Db) given 
U G SO(Da) in the other Franson loop. Since SO(D) is the endomorphism group of the 
D-dimensional redit state space we may consider this restriction of the maximization 
procedure when the state space is restricted to a two-redit subspace. For two redits the 
state p naturally decomposes as 



XDl- 



P 



[5 ko (D A -2)+2)[5 l0 (D B -2)+2] 



EE(* 

k=0 1=0 



sym 
kl 



(27) 



where o fe " 



is nonzero only when xt an d xf are both symmetric, and S\ 



tanti—sym 
'kl 



IS 



nonzero only when \k an d xf are both antisymmetric. The generators of the subgroup 
SO(D) C SU(D) are the antisymmetric generators of SU(D). As a special case we 
can consider a general mixed two-rebit state, where the only pair of antisymmetric 
generators spanning the state space is Cgxr-f . Since SO (2) operators can be expanded 
in a basis consisting of only 1 and <r 2 , the intensity is 

jAB 



~ + \u Q V - -V 2 U 2 M 22 , 



(28) 



where U , V , U 2 , V 2 e R and M 22 = C R 
The SO (2) operator that maximizes the intensity is 

V = j (U 1 B - iC R U 2 a B ) , 



a B p) is the rebit concurrence 



(29) 
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where A = v/j Uq \ 2 \ U2 | 2 . We note that for product two-rebit states V can only be 
1 or — 1. 

4. Correlation induced non-Abelian quantum holonomy 

In this section we use the parallelity condition introduced in section 13.11 to define 
a procedure for parallel transport of a bipartite quantum state. We consider the 
infinitesimal limit to find a connection form corresponding to this parallel transort. 
By definition the output of the long arms is parallel with the output of the short arms 
when V is chosen such as to maximize the coincidence intensity. Now we choose to view 
the output of the long arms as the parallel transported version of the output of the 
short arms. By using that output state as input for another Franson setup, where in a 
similar way a new output state is created, we can parallel transport the state through 
an arbitrary number of steps. 

To see how this works we let p^ be the input state of the interferometer. The 
output state of the long arms is p (1) = U w ® V^p^U^ <g) \/ (1)t , where G U(D A ) 
and V W G U{D B ) are unitary operators that has been applied such as to implement 
parallelity. In the second step, we use p^ 1 ' as the input in a new Franson interferometer, 
where a new unitary U^> is chosen and a new is found to create an output of the 
long arms p*- 2 -* that is parallel to p^\ 

The parallel transport is performed by iterating the intensity maximizing procedure 
in this way as illustrated in Fig. 2. In the nth step, a G U(Da) is chosen, and 
thereafter a G U(D B ) is found that maximizes the intensity. After has been 
found the input state for the next step is taken to be p^ = (U^ (&V^)p( n ~ l \U^ ® 
V'Mt). 

The coincidence intensity in the (n + l)th step is now 



J (n+1) = I + ^ReTr [(f> (n+1) ®y (n+1) )(^(n)®yH)p (0) (?7(n) t ®y(n) t ) 



2 2 



(30) 



where the cumulated unitary operations that are applied to the original input state 



(0) 



at the beginning of the (n + l)th step are U{n) = U^U^ n X K . . and V{n) 



P 

y{n)y{n-l) _ _ _ y(l)_ 

From this we can define a holonomy group Hols based on a particular state 
p with corresponding Stokes tensor S as the set of unitary operators V(n) = 
y(n)y(n-l) _ _ y(i) e U(D B ) that can result from the above parallel transport 
prescription given all sequences of G U(Da) such that U (n) = (ji^U^- 1 ) . . . (jw = 
1 for any n. From the discussion on product states at the end of section 13.11 follows that 
the holonomy group for product states is always Abelian, and only correlated states can 
induce a non-Abelian holonomy group. 
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Figure 2. Iterative procedure in first, second and (n + l)th step. 
For any set of unitaries {U^ G U(Da)} given by 

D\-l 

U(n )=u (n )iA+ J2 Utht (31) 



k=l 



we find G U(Db) that maximizes the intensity as 



j,fc=0 1=0 



(32) 



where = Tr 



(xf®x£)(U(n - \)®V{n - l))p°(U(n - l)i®V(n - l) f ) . 

To define a connection we need to consider the limit when U and V are 
infinitesimally close to unity. To find this limit we revisit the maximization problem 

with a different parametrization U = e l ^>=° and V = e*^=° . The intensity 

is then 

I AB = \ + ^Re Tr(e* ^A~\xf ® e * E 5~^^p ) , (33) 

where flj- and 0-,- are real numbers. In this representation unitarity is explicit and no 
constraints are necessary. Differentiating I AB with respect to the parameters 4>j an d 
setting each derivative to zero, we find 

ReTr(e iE ^ l0fc ^ ® iS^e^^'^po) =0, < j < D| - 1. 
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(34) 

Since we are only interested in finding the connection form we expand these equations 
to linear order in 9j and (f>j, to obtain 



D\-l 



D|-l 



ReTr \l A +i 9 k 



<PixF I Po 



k=0 



1=0 



- { ^2<PkS k5jo + (j)oS 0j + 



k=0 



D, 



-^+1-^)6 



>j0 



<i>j / - y^ <i>kdjkiSoi 



k,i=i 



D\-l 



y^ ^'S'fcj — o, 



< j < D% - 1. 



(35) 



fc=0 



In the infinitesimal limit we introduce the notation dUU^ = i(52f=o * dOjxf) 



and likewise dVV^ = i(X]^o 1 d(j)jxf)- Although performing infinitesimal unitary 
operations is clearly an idealization we may still consider this limit where the sequences 
of unitaries {A^ } = 1 A + dUU\t)} and {AJ } = 1 B + dVV\t)} in the parallel transport 
are indexed by a continuous variable t. 

The relation between {A^} and {Ay }, for each t, is given by 



D|-l 



E B jk (dVV%(t) = -J2 (dUU%(t)S kj , 



(36) 



fc=0 



k=0 



where Bj k are the elements of the symmetric matrix B, given as 



B, k = ReTr(xfx!p B ) 



Dr 



5 kj + ( 1 - 



D, 



Dl-l 



+ E S ° l [(^'o^fc + SkoSij) + d jk i] 



i=i 



Therefore, provided B is invertible, we find 



D\-ID%-1 



= i B - E E (duu%s M B-]x*. 

k=0 l,m=0 

We can now identify 



(37) 



(38) 



Dl-lDl-l 



A ® = ~ E E (d^)fc5„S^ (39) 

fc=0 l,m=0 

as the operator-valued anti-Hermitean connection one-form. Note that A is a linear 
function of the Stokes matrix S. If we decompose the density operator as p = Yl^P^ I 
VvK^V I we can express the connection form as 



D%-1 

A(t)=iJ2p»J2 ReTr ( 1A ® i xf \d^,{t))^,{t)\)Bj l X f. (40) 

H j,l=0 
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Under a change of gauge | ipk) ~> 1 A ®G \ ipk), corresponding to a unitary transformation 
on the second subsystem, the connection transforms as 

D|-l 



M j,l,m,n=0 
D%-1 



z^^ReTr[l A ®Gt^B ( G \d%) + dG l^^R^R 

j,l,m,n=0 



T ~B 
mnA.n 



M h,j,l,m,n=0 
D 2 -1 



+*E*V E ReTr[l A ®^fGtrfG | ^)(^ \]R T hj R 3l B^R T mn x* 

fi h,j,l,m,n=0 

= GdG ] +GAG\ 



(41) 



where i? ifc = feo(gfl 1 - 2 ) +2 ] Tr (^xf G*xf ), and we have used that B jk = ReTr(xfxfp B ) 
which transforms as B — > RBR T . Thus A transforms as a proper gauge potential. 

For a given path 7 in U(Na) given by £/(t), the parallel transport gives us a path 
V(t) in £/(JV B ) 



V(t) 



exp 



A(s)ds 



(42) 



where P denotes path ordering. The holonomy for a closed path in U(Na) is thus given 
by such an integral and is dependent on the Stokes matrix via the connection form in 
equation (1391) . 



5. Relation to Levay parallel transport for SU(2) x SU(2) 

The pure two-qubit states can be represented as quaternionic qubit states [23], [21] using 
the structure of the second Hopf-fibration. Within this representation one can construct 
the quaternionic analogue of the Pancharatnam geometric phase, as has been done by 
Levay [21]. We review this quaternionic representation and show that when the state 
evolution is generated by local SU (2) operators, the Levay geometric phase is contained 
in our construction. 

In the quaternionic representation a pure two-qubit state 

I ^) = a 1 00) + p 1 01) + 7 1 10) + 6\ 11), (43) 
where a, ft, 7, 5 G C, is associated with a quaternionic qubit state 



[a 



m io)+( 7 +fi) ii>. 



(44) 
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Here the quaternionic state | 0) is formed by identifying | 00) and | 01) in such a 
way that the complex coefficients of these states are contained in the new quaterionic 
coefficient, and similarly for | 1). The standard quaternion basis elements i, j and k, 
satisfy i 2 = j 2 = k 2 = — 1 and ij = — ji = k. A local SU(2) unitary U acting on the 
first qubit is represented by a operator U acting from the left on | and a local SU{2) 
unitary V acting on the second qubit is represented by a unit quaternion q v acting from 
the right 

|*) = U|* )gy, (45) 

where U = Uol + iY^=i^j^j with i and <ti,<t 2 ,<X3 the standard unit and Pauli 
operators acting on the quaternionic qubit Hilbert space. The unit quaternion q v E 
Sp(l) corresponds to V according to 

3 



VjC' B 



t 

q v = V + V 3 i - V 2 j + V x k, (46) 

where a B ,a B ,a B are the standard Pauli operators acting on the state space of the 
second qubit. 

The inner product of two quaternionic states \ = p± | 0) + q± | 1) and 
I*) = P2 | 0) + q 2 1 1) is 

(*|$) = qlq 2 +p* 1 p 2l (47) 

where * is the quaternionic conjugation operation defined by (a + bi + cj + dk)* = 
(a — bi — cj — dk), a, b,c,d E M. 

The quaternionic transition amplitude between two states related by a local SU(2) 
operation U on the first qubit can be expressed in terms of transition amplitudes in the 
ordinary complex representation as 

3 3 3 

i u i *) = u + J2 u j m ^ - E u s M & + E U ^ M ^ ( 48 ) 
j=i j=\ j=\ 

where | ip) and | ^) are the complex and quaternionic representations of the same 
state. From this we can consider the formal analogue of the Mach-Zehnder interference 
intensity 

I= l -+ l -Re{(*\t\mv) 

= \ + \u V - ifpjVkMjk, (49) 
j,k=i 

where the phase factor qv represents a local SU(2) operation on the second qubit. To 
compare this intensity to the Franson interference intensity for the case U E SU(2) 
and V E SU(2), we consider equation (124"]) when Da = 2 and use the parameterization 
U = Uq1 a + i J2 3 j=i Ujof, where Uq, U\, U2, U3 are real numbers. We then find that the 
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quaternionic Mach-Zehnder interference intensity J, is identical to the Franson intensity 
I AB , which demonstrates the correspondence between the quantum- mechanical Franson 
setup and the quaternionic quantum-mechanics Mach-Zehnder setup as shown in figure 

El 



U 



a — L \ o / vh) 



V 



u 



4) 



Figure 3. Correspondence between Franson setup and quaternionic QM Mach- 
Zehnder setup 

This quaternionic representation of the two-qubit states can be considered a 
generalization of the ordinary complex spinor representation of single qubits. In the 
ordinary single qubit representation, the Hilbert space is the space of normalized 
complex spinors S 3 and the projective Hilbert space is the space of spinors modulo 
a phase factor S 3 /U(l) = S 3 /S 1 = S 2 = CP 1 , the complex projective space of 
complex dimension one. In this quaternionic representation of two-qubit states, the 
Hilbert space is the space of normalized quaternionic spinors S 7 . The quaternionic 
projective Hilbert space is the space of quaternionic spinors modulo a unit quaternion 
acting from the left, representing, as mentioned above, a SU(2) = S 3 rotation of 
the second qubit. This implies that the quaternionic projective Hilbert space is 
S 7 /S 3 = S A = HP 1 , the quaternionic projective space of quaternionic dimension one, or 
real dimension 4. While the ordinary single qubit representation corresponds to the first 
Hopf fibration S 1 "-^ S 3 "-^ S 2 , the quaternionic two qubit representation corresponds to 
the second Hopf fibration S 3 "-^ S 7 "-^ S A . This quaternionic two-qubit representation 
is much similar to a qubit in quaternionic quantum mechanics [31] except that in this 
representation the absolute quaternionic phase corresponds to a rotation of the second 
qubit and thus is a measurable quantity. 

Within this representation of pure two qubit states, Levay [24J studied the 
quaternionic analogue of the Pancharatnam parallel transport. Levay's parallelity 
condition and related parallel transport are defined such that two quaternionic states 
are parallel if their inner product is a real and positive number. This condition is in 
concordance with the Mach-Zehnder analogue picture since the maximal intensity is 
achieved when the unit quaternion phase factor q\ is such that | U \ *&)qv is real and 
positive. We next consider a family of spinors that are stepwise unitarily evolved from 
the same initial state | \I/o) and demand that in each step the initial and final states are 
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parallel. The transition amplitude between two consecutive spinors in this family is 
| * n > = (* n _! | UW | *„-i)g£° = A, AG R+, 

(50) 

where | Vf^-i) = U(n — 1) | & )qv(n — 1), and we use the notation U(n — 1) = 
I j(n~i) I j(n-2) _ _ _ and q v (n - 1) = ffiffl . . . q ( y~ l) . Thus 

(*„_i|U< n > |* n _ 1 ) = Ag? ) *, (51) 
or by using equations (jlBjl and 



Hence 



£/ (n) + $>j n) M j3 i -£>j n) M j2 j +^[/f ) M il k = 
j=i j=i j=i 

= X(V Q {n) - V^h + V 2 {n) j - V{ n) k). (52) 

T/M _ ^ rr( n ) 

Vf ] = - \Y.U^M kr (53) 

k=l 

The parameter A must be chosen to normalize V^ n \ hence A 2 = t^o 2 +X}i(Sj U^Mji) 2 . 
If we compare this to equation (|25|) . when Da = 2 and again use the parameterization 
U = Uol A + iY^=i^j^fj where Uo,Ui,U2,U3 are real numbers, we find that this 
parallelity condition is the same as that in equation (1531) . 

To find the Levay connection, we consider the infinitesimal limit in which the 
parallel transport condition reads 

(d* | *) = 0, (54) 

where | \&) is the instantaneous state. If we only allow changes generated by the local 
unitaries U and qy we have 

(d* | *) = dq v *q v + (d& | UdU f | (55) 

Imposing the parallel transport condition (d^f \ = 0, and using equations (jlBj) and 
18]) we find this to be equivalent to 



(dVV^j = -J2( d UU r ) k M kj . 



(56) 



k=i 



To compare this expression with the connection in our construction we consider equation 
(I36p for Da = Db = 2, and note that for SU(2) all symmetric structure constants djki 
are zero. This gives the B matrix the following form 

/ 1 iSm Stv) Srra \ 



B 



'01 °Q2 ^03 

S 10 1 

S 20 1 

S 30 1 



(57) 
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By taking into account that (dUW) = and (dVV^) = for 5*7(2), we see that only 
the correlation matrix M will be relevant to the relation between dUU^ and dVV . Since 
Bjk = Sjk for j, k 7^ it immediately follows that equation ( 1361) reduces to equation ( 1561) . 

6. Conclusion 

We have constructed a parallel transport procedure in the same spirit as that of 
Pancharatnam [HI IB], in the sense that it defines parallelity with reference to 
maximization of an interferometric quantity. The interferometric quantity chosen in 
this case is the coincidence intensity of a Franson type interferometer. The phase is 
taken to be the local degrees of freedom of one of the subsystems. 

Given two different two-partite states related by local unitary evolution of one 
of the subsystems, the unitary operation that needs to be applied to the other 
subsystem to achieve parallelity, depends on the correlation present in the full bipartite 
system. Generally phase unitaries that correspond to different parallel transports do not 
commute, however, when the system is uncorrelated, only Abelian phase factors need 
to be applied. Thus the holonomy group related to the parallel transport condition is 
Abelian if the bipartite state is uncorrelated, and a non-Abelian holonomy group can 
be said to be correlation induced. 

The procedure is defined for arbitrary bipartite systems, pure as well as mixed. 
In the infintesimal limit of the parallel transport, the connection form can be found 
as a closed expression for arbitrary dimension. On the other hand, finding a closed 
expression for parallel transport when the steps are finite appears to be non-trivial in 
the general U(D) case. The procedure can be restricted to subgroups of the full unitary 
groups. In the pure two-qubit case when only SU(2) operators are considered it has 
been shown that this procedure is related to Levay parallel transport |24j. Therefore 
our construction opens up for experimental tests of the Levay geometric phase in the 
special case of local SU(2) evolutions. 
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